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ABSTRACT 

The long-time large-distance behaviour of free decaying two dimensional turbulence is studied. 
Stochastic solutions of the Navier-Stokes equation are explicitly shown to follow renormalisation 
group trajectories. It is proven that solutions of the Navier-Stokes equation asymptotically converge 
to fixed points which are conformal field theories. A particular fixed point is given by the free 
Gaussian field with a charge at infinity. The stream function is identified with a vertex operator. 
It happens that this solution also admits constant n-enstrophy fluxes in the asymptotic regime, 
therefore fulfilling all requirements to represent an asymptotic state of two-dimensional turbulence. 
The renormalisation basin of attraction of this fixed point consists of a charged Coulomb gas. This 
Coulomb gas gives an effective description of turbulence. 
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I Introduction 

Two-dimensional turbulence is a simplified version of three-dimensional turbulence. It can 
nonetheless give qualitative information on phenomena such as atmospheric turbulence. Recently, 
many field theorists have been attracted by 2d turbulence. This new interest has been aroused 
by Polyakov's seminal' 1 ! paper. For the first time, it seems that this long-standing problem can 
be tackled using field theoretic methods' 2 !. Definite predictions are obtained. One can hope to 
compare them with numerical simulations' 3 ]. 

Most analyses of two-dimensional turbulence aim at computing the two-dimensional analogue 
of the celebrated Kolmogorov spectrum' 4 ] in the long wavelength limit . An early analysis by 
Kraichnan' 5 ! showed that the energy spectrum should decrease as A; -3 , due to the transfer of 
vorticity from large to small scales where dissipation takes place. Nevertheless, this prediction 
suffers from large infrared divergences and is corrected by logarithms at one loop order. On the 
other hand, such a correction does not assure convergence at two loop order where double logarithm 
infrared divergences appear. Alternative predictions have been given by Saffman' 6 ! and Moffat' 7 ! 
with respective spectra k~ 4 and k~~ . This variety of results shows that the situation is far from 
being settled. From a field theoretical point of view, these approaches have severe drawbacks, they 
all assume that dimensional analysis is exact and neglect possible non-trivial anomalous dimensions 
appearing in short distance expansions of field products. Unfortunately, no precise rationale can be 
extracted from numerical simulations. It happens that simulations starting from gaussian initial 
conditions show an increase of the exponent from -4 to -3 for short times. Then coherent vortices 
emerge. Moreover, a very long time simulation' 8 ! shows that the fluid motion seems to converge 
to an ensemble of point vortices' 9 !. The fate of the spectrum exponent as time goes to infinity is 
certainly not clear. 

The Polyakov approach deals with steady states of two-dimensional turbulence. This amounts 
to imposing boundary conditions and an external pumping force. Turbulence is then only obtained 
in the inviscid limit where the viscosity goes to zero. Interpreting the Navier-stokes equation as an 
equation for correlation functions, Polyakov shows that scale- invariant field theories can represent 
steady states if a certain inequality is satisfied by operator dimensions. This boils down to imposing 
a vanishing non-linear term in the Navier-Stokes equation. Another constraint is obtained by 
forcing the enstrophy flux to be a constant. A simplified treatment of the resulting Diophantine 
equation leads to a solution whose spectrum scales as k~ 25 ^ 7 . Unfortunately, such a steady state is 
non-unitary and thus suffers from large infrared problems. The Diophantine equation can be shown 
to yield thousands of possible non-unitary minimal conformal field theories' 10 !. As emphasised in 
ref.fll], these solutions are hardly physical. In addition to infrared divergences, they fail to provide 
a non-zero enstrophy flux. Moreover, they only take into account the enstrophy flux and discard 
fluxes of higher conserved quantities. It has been recently shown ' 12 1 that a conformal description 
of two dimensional steady states is rather unlikely. 

In this paper, we strive to tackle a less ambitious problem, i.e free decaying two dimensional 
turbulence. To do so, we follow closely the Polyakov approach but try to remedy some of its 
drawbacks by imposing that fluxes of all conserved quantities are constant. As time increases, 
viscosity effects become less important and solutions of the Navier-Stokes equation reach a quasi- 
equilibrium containing free coherent vortices. For such a state, it seems reasonable to expect that 
fluxes as well as the energy are zero. The object of this article is to show that these conclusions 
can be derived from an effective field theory describing the long-time regime of two-dimensional 
turbulence. 

The paper is arranged as follows. In a first part, some comments on the nature of turbulence 
are proposed. The second part deals with solutions of Navier-Stokes equations. In the third part, we 
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impose that fluxes of conserved quantities are constant. Then, the long time behaviour of conserved 
charges is analysed. In the last section, we propose an effective theory for two dimensional free 
decaying turbulence. 



II A Field Theoretical View of Turbulence 



Turbulence is often described as a random behaviour reached by fluids under certain conditions 
such as high speed and low viscosity. Fluid dynamics describe any fluid motion in terms of the 
velocity field v. This field is a macroscopic field: it is the average of molecular velocities over a 
macroscopic size a ( a is rather small and can be thought to be of the order of a micron). Below 
this scale, the Hydrodynamic approximation is not valid. There is another intrinsic scale Z; it is 
set by viscosity effects and represents the scale below which viscosity smoothes variations of the 
velocity field. Above this scale, the velocity is rapidly varying. A sufficiently low viscosity is a 
phenomenological condition for the existence of a turbulent regime. In that case, the viscous length 
I becomes of the order of the macroscopic length a. We suppose that this criterion is fulfilled, so 
that there is a single fundamental length scale a. 

The equation governing the time evolution of the velocity field is a non-linear equation, 

(d t + v.d)v = uAv ~Vp + F (1) 

known as the Navier-Stokes equation where p is the pressure field and F the external pumping 
force. Such an equation is simpler in two dimensions as it only involves a pseudo-scalar field 
instead of a vector field. This allows strong results about classical solutions' 12 ] to be proven. It 
can be shown that smooth (C°°) initial conditions remain smooth. However, this is not the sort of 
statement we are interested in. Indeed, we would like to understand the Navier Stokes equation 
when the velocity field is far from being smooth. 

As a matter of fact, the velocity field is a random field. This means that each field configu- 
ration is chosen at random according to a probability law. We will restrict ourselves to probability 
laws which can be expressed as statistical field theories. Finding the probability law amounts to 
finding the action S t giving the Boltzmann weight at every time t. The action S t is an effective 
action valid for distances larger than the macroscopic length a. It is a local action. Indeed, the 
range of the interaction between two fluid elements (by this we mean a macroscopic neighbourhood 
of a point) is given by the viscosity length I. Below this scale, the velocity field hardly varies due 
to viscosity effects; on the contrary, for larger scales, two fluid elements are not directly in contact, 
and their velocity fields can be almost uncorrelated. But for turbulent fluids, the viscous length 
is of the order of the macroscopic length. The length a plays the role of a physical Ultra- Violet 
cut-off, so that the interaction length is of the order of the UV cut-off. As all fields will be smeared 
out (regularised) over a distance a, this implies that the Lagrangian composing the action is in 
fact local, i.e. can be expressed as a function of regularised fields and their derivatives. This does 
not prevent the existence of any long-range order. As for critical phenomena, long-range order can 
occur due to collective effects even though interactions are short-ranged. 

Two highly different issues could be tackled. On the one hand, the Ultra- Violet problem 
amounts to solving the Navier-Stokes equation with an external pumping force and given boundary 
conditions. Most of the time, a large box and a low frequency random pumping force are specified. 
Turbulence only appears in the inviscid limit where v goes to zero. The velocity field is random 
because of the random external force. Notice that in field theoretic parlance, this is tantamount 
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to studying the UV behaviour of the action S t , keeping an Infra-Red cut off fixed. Scale invariant 
solutions are then time-independent UV fixed points in the renormalisation group sense. They 
represent stationary turbulence states. On the other hand, the Infra-Red problem deals with 
infinite space and Dirichlet boundary conditions at infinity. No external pumping is required but 
random initial conditions are imposed. As shown in section III, the viscosity v is finite and does 
not need to be scaled down to zero. It becomes negligible in the long time regime. This is a free 
decaying turbulence situation where scale invariance is gradually restored by time evolution. It is 
quite analogous to a self-criticality regime. This issue has to be faced as an Infra-Red problem, 
where long-time, large-scale limits are taken, keeping the UV cut-off a fixed. Asymptotic solutions 
are then IR fixed points. The relation between IR and UV fixed points seems to be a difficult issue. 

It is convenient to move from the velocity picture to the vorticity picture. In 2d, the continuity 
equation reads 

Vv = 0, (2) 
valid for an incompressible fluid. This entails that v is a curl: 

Va = e a /3dpip (3) 

where ip is the stream function. The field -0 is a pseudo-scalar field. In the plane, it is always 
possible to construct pseudo-scalar fields from scalar fields, they are just two-forms ipdz A dz. It is 
also useful to introduce the vorticity which is the curl of v: 

u = e a/ 3d a vp. (4) 

Notice that ei 2 = — e 2 i = 1 is the antisymmetric tensor, i.e. u> = —Aij). This is also a pseudo-scalar 
field. This yields a different form of the Navier-Stokes equation: 

d t u + e a pd a ipd[3Aip = vAlu (5) 

This equation for the vorticity field will be the main ingredient of our analysis. 

Before further delving into the properties of the vorticity field, we will recall a few facts about 
certain two-dimensional field theories. The velocity field obeys the Navier-Stokes equation. It is 
not clear what meaning can be given to the product of fields appearing in the non-linear term. 
The following will provide an answer. Fields will be chosen to be tempered distributions G S' 
(continuous linear forms on the space of rapidly decreasing functions S). A field theory such as 
the IR fixed point is a measure on the space of tempered distributions' 13 ^ To regularise the short 
distance behaviour of fields, introduce a smoothing function p e : 



I 



d 2 x p(x) 



f \ 1 f x \ ( 7 ) 

Pe{x) = -^P(-) 

P > 0. 



This function is an approximate delta function. It allows to transform the very rough variations of 
tempered distributions into perfectly well behaved functions. To do so, introduce the regularised 
field 

4> e (x) = J d 2 yip(y)p e (y-x) (8) 
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This is now a smooth (C°°) function. We will interpret macroscopic fields as regularised fields 
where e = a the physical UV cut-off. This allows to see the Navier-Stokes equation and its vorticity 
counterpart as equations for regularised fields. The non-linear term is then a simple product of 
smooth functions. Obviously care will be needed in removing the cut-off. 

The Navier-Stokes equation is now expressed in terms of regularised fields. Physically, this is 
all the more natural as the velocity field is the average of molecular velocities. This gives: 

d t uj a + e a pd a u a df3Aip a = uAu a . (9) 

The initial time field theory is specified down to the cut-off a. It is the initial state of the Cauchy 
problem for the evolution equation. It is specified by an initial effective field theory representing 
the initial probability law of the stream function. At any given time t, the effective field theory 
describing the behaviour of ip a is obtained from (9). It represents the probability law describing 
the fluctuations of ip a at time t. 

Suppose that the effects of viscosity can be neglected in the asymptotic regime. Then the 
inviscid fluid possesses an infinite number of conserved charges. The first one is obvious, it is the 
energy per unit mass. Another set of conserved quantities stems from the pseudo-scalar behaviour 
of the vorticity. Introduce the generalised n-enstrophy: 

H n , a = J d 2 x u2(x). (10) 

Only the even enstrophies are non zero. Now compute 

= ne a/3 J d 2 x daipadpAipa LOa~ x (x). (11) 

Integrate by parts using the antisymmetry of e a fs to obtain: 

z a p J d 2 x uj™d a dpAuja = 0. (12) 



dH na 

e c 



dt 

We will need to integrate by parts very frequently. So from now on, we restrict fields to be zero 
at infinity, ie regularisations tp a , u a ... vanish at infinity. The second moment is simply called the 
Enstrophy. In section III, we will justify this calculation directly from the viscous Navier-Stokes 
equation. It will be shown that this result is indeed true in the long-time regime, even in the 
presence of viscosity. Such conserved charges give rise to conserved currents in Fourier space. 
These currents represent the flux of n-enstrophies from scale q to q + dq. In the usual Kolmogorov 
approach, the energy flux (called the energy transfer rate e) is supposed to be finite and constant. 
Physically, these conditions imply that every point in the fluid plays the role of a local source of 
n-enstrophy and energy. In the free decaying case, we expect these fluxes to be zero if a steady 
state exists. In the following we will call turbulence any critical solution of Navier-Stokes equations 
such that the associated fluxes are finite and constant (possibly zero) . 
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Ill General Asymptotic Properties Of Decaying Turbulence 



a) The Neighbourhood of Renormalisation Group Fixed Points 

We will suppose that turbulence is specified by a massless field theory. This theory is a fixed 
point of the evolution equations in the long time-large scale regime. From given initial conditions, 
it is not clear what the set of fixed points can be. It could be a single theory, and in that case any 
initial condition would finally be attracted by this fixed point. It could also be a more complicated 
set such manifold. 

On physical grounds, 2d turbulence is likely to be represented by a unique fixed point. Indeed, 
2d turbulence is isotropic. Then, universality suggests that there must be a unique universality class 
describing 2d turbulence in the infrared regime. This is certainly true if the notion of universality 
can be extended to self-organised criticality. In the following, we will try to first analyse the features 
of a fixed point. Thus, we will show that any conformal theory satisfying a specified inequality 
can be a fixed point. Fixed points are actually asymptotic solutions of Navier-Stokes if they are 
stable fixed points, i.e. if nearby deformations of a fixed point converge in the renormalisation 
group sense. 

First observe that the viscosity defines a length scale and a time scale, i.e. by engineering 
dimensional analysis: 



a 2 



(13) 



The scale a is the UV cut-off whereas typical times are measured in units of r. We want to 
investigate the large distance r 3> a long time t>r regime. Let us first define a fixed point as a 
continuum conformal field theory satisfying the Euler equation: 

lim [ f + e a pd a i>«.dpAiPa] =0. (14) 

Notice that such an equation involves the product at shorter and shorter distances of two regularised 
fields. In order to further analyse the Euler equation, we need to evaluate short distance expansions. 
In a nutshell, the expansion 

ii(y)4>(z) = J2c a (y,z)O a (z) (15) 

a 

is valid in any correlation function. The coefficient functions are in fact distributions which scale 
as v = \y — z\ — > as: 

C a (y,z) ~ i/toc.-**-**- (16) 

Moreover, the most singular term on the right-hand side can be retained as the most significant 
piece. The Euler equation involves a non-linear term whose continuum limit is far from being 
obvious. Nevertheless, it can be computed using the operator product of two tp. Indeed, the 
non-linear term reads (we put e = f ) 

e a i3d a ^e{x)dp/\^ e {x) = j d 2 ud 2 ve a/3 d a p e (u - x)dpAp e (v - x)ip(u)ip(v). (17) 

In the continuum limit e — ► 0, it is safe to retain the dominant term in the operator product 
expansion: 

P d 2 ud 2 ve af 3p € (u - x))p £ (v - x)d^A v J2C 1 (u,v)^(v), (18) 
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where ^2 are the fields appearing in the short distance expansion. The only way of getting a 
non-zero term is to act on the possibly non-symmetric term under a — > (3. One of the coefficient 
functions must satisfy C(ui,ii2,Vi,V2) = —C(v 2 ,vi,U2,ui). This can only be obtained from de- 
scendent terms, where the leading operator is constructed from a pseudo-scalar combination of 
descendent fields L— n i....L- nm ip. The Virasoro generators L n = J dzz n+1 T(z) are the gener- 
ators of the 2-d conformal algebra and T is the energy momentum tensor. Truncating the series 
yields: 

J d 2 vd 2 up e (v) Pe (u-x)e a(3 d v a d%A u C 1 (u + v,u)Mu) (19) 
After changing variable v — > ev, and using the fact that the cut-off is small: 

J d'vp^e^dld^C^u + ev.u) ~ ^r 2 ^ -4 J cPve a p^A u C^u + v,u). (20) 

We will denote the constant by C, neglecting its dependence on u: 

e^d^dpA^ ~ e d ^~ 2d ^ 4 C^l e (x). (21) 



From now on, all fields will be normalised in such a way that finite constants like C are equal to 
one. It is now necessary to identify the leading operator. Notice that due to the term e a p the 
result is a pseudo-scalar. Futhermore it has to be antisymmetric under u\ — > u 2 . Using complex 
coordinates: 

^2(z,z) = -^(z,z) 



Under the first transformation the energy momentum tensor transforms as: 

(23) 



T -» -T 
T — > —T 



This yields the transformations of Virasoro generators L n — > i n L n . The lowest order anti-invariant 
combination is: 

iP 2 = {L_ 2 L\ - L_ 2 L 2 -M2. (24) 
This immediately gives the leading behaviour in (20) as 

d^ = d^ 2 + 4. (25) 

This calculation can be extended by induction to higher order short distance expansions. After 
this lengthy computation, we obtain the conditions for the existence of a fixed point. Suppose that 
the fixed point theory is time invariant. Then the Euler equation is satisfied if the non-linear term 
vanishes altogether in the continuum limit. This implies that' 1 ': 

di, 2 > 2d i , (26) 

This is the Polyakov inequality. Hence, all conformal field theories satisfying this equation are 
solutions of the Euler equation. However, they may not be valid representation of 2d turbulence. 
To select the expected unique fixed point modelling 2d turbulence, we will have to impose that 
fluxes of n-enstrophies are constant in the infra-red regime. 
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We are now in position to prove that fixed points are asymptotic solutions of the Navier- 
Stokes equation. Suppose that, after a transient period of order O(r), the field theory obtained 
solving the Cauchy problem from a given initial condition falls within the renormalisation group 
neighbourhood of a fixed point (we will call this neighbourhood the basin of attraction of the fixed 
point). This means that this field theory converges in the large distance limit towards the fixed 
point. Let us prove that such renormalisation group trajectories coincide with the Navier-Stokes 
evolution. Suppose that large distance long time trajectories are renormalisation group trajectories. 
To implement this idea, rescale x — > Ax and t — > \ T t where T is a yet unknown constant. Then 
define a new field ip as the regularised field following the renormalisation trajectories of the solution 
ijj at time t. This means that the regularised field V> a (Ax, X T t) obeys the following: 

M\x,\ T t) = \- d ^ x - la (x,t). (27) 

The field tp is the effective field obtained after integrating the action St over distances in the range 
[A _1 a, a], whereas the field ip\-i a is the regularised field defined with the action S t - Notice that 
this ansatz ^ is a solution of Navier-Stokes equation provided that, when plugging ip in (9), one 
ends up with the original the Navier-Stokes equation for ip at time t and cut-off j . Obviously, such 
an equation should have a modified viscosity v — > a2 -t as the typical time r is also rescaled to w • 
Putting this ansatz back in Navier-Stokes equation after adjusting the time rescaling t — > X 2+d ^t, 
one obtains 

d t u\-i a + e a pd a ^ x -i a dp^ x -i a = X~ 2+T uAu x -i a . (28) 

This is indeed the expected result. This proves that -0 is a solution. Notice that scale invariance is 
asymptotically achieved if and only if the dimension is negative. Moreover T has to be positive 
in order to ensure that time flows towards infinity. This implies that —2 < < 0. Hence, we 
retrieve Polyakov's idea that turbulence requires non-unitary theories. This entails an extreme 
boundary condition sensitivity. We will end up considering charges at infinity which in a sense 
epitomise the non-trivial infrared properties of turbulence. Notice that viscosity explicitly breaks 
scale invariance. The effective action in the basin of attraction of the fixed point will reflect this 
fact, i.e. it will be the sum of the fixed point action and irrelevant operators flowing to zero in 
the infrared regime. To conclude, the field ip a \-i tends to the continuum limit field -0 whereas the 
solution tp converges to the fixed point along renormalisation group trajectories. This proves that 
fixed points are asymptotic solutions of Navier-Stokes equation 1 . Amongst all these fixed points, 
we will explicitly find an example admitting finite fluxes in the infra-red regime. 

c) The Effects of Viscosity on iV-Enstrophies 

As viscosity is introduced, it is no longer true that ra-enstrophies are conserved. Nevertheless, 
they are slowly varying. In order to assess the effect of viscosity, let us compute the ra-enstrophies 
and their time derivatives in the asymptotic regime. By (10), 

H nta (X T t)= [ d 2 x^(x,A T t), (29) 



1 This result is also valid in three dimensions. In that case, solutions of the Navier-Stokes 
equation still follow renormalisation group trajectories. In the infrared limit, solutions converge 
to non-unitary conformal fixed points defined by a straightforward generalisation of (14). On the 
contrary, the Polyakov inequality seems to require a non-trivial extension. 
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Change x — > Ax and use the scaling properties of u> to obtain 

^n,a(A T t) = A 2 " 2 "-^ | d 2 X l^-i (x,t). (30) 

The integrand tends to be the product of n fields at the same point, so it is necessary to evaluate 
higher order short distance expansions. These products are defined by induction. They can be 
written in a compact way using conformal families: 

[#Vn-l] = [Vg + .... (31) 

where tp n = Dom(-0VVi-i) is the leading part of the short distance expansion ipip n _i. This entails 
that inductively: 

w"-i„0M) ~ (Aa- 1 )-«**»+" d *+ 2 "V niA -i (x,t). (32) 

Hence the n-enstrophy: 

H n , a {\ T t) = A 2-^ no ^-n^-2„y d 2 x ^ n A _ la ( X)t ). (33) 

./V-Enstrophies are simply integrals of higher order fields. 
Similarly, the time derivative reads 

dHn 'dt XTt) = n ^~ 2n ~ n ^ J d2x ^Au; x -i a (x,t)u;^ a (x,t) (34) 

where the rescaled version (28) of Navier-Stokes equation is used. The right-hand side is always 
negative as integration by parts yields — (n — l)u j d 2 x (9cj A -i a ) 2 oj™r 2 a . Finally, 



dH n , a (\ T t) 



-n^A 4+ ^- d ^a- 2(n+1) - n ^ + ^« J d 2 x ip n ,x-i a (x,t). (35) 



dt 

Using the scaling properties of the fields, this yields in the infrared regime, 

= ~ n ^At) (36) 

where t is now the asymptotic physical time. Notice that the transient time ^ appears. The 
re-enstrophies go to zero in the asymptotic regime. Similarly, one gets 

E a (t) = £ a ( r )exp(-2— ) (37) 

r 

in the asymptotic regime. Notice that the energy is exponentially damped to zero unless there is 
no viscosity. 

So we have found that in the infrared regime solutions of Navier-Stokes equation have vanishing 
conserved quantities. This does not imply however that solutions are trivial. 



c) The Energy Spectrum 
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In the same vein, one can obtain an expression for the scaling exponent of the energy spectrum. 
The spectrum is nothing but the Fourier transform of the velocity connected 2-point function: 

E(k) = 2irk J d 2 xexp(— i2irk.x) < v a ,a(0)v a , a (x) > (38) 

such that 

E a = j dkE(k) (39) 

The scaling behaviour of E(k) in the inertial range can be deduced as follows: 

< d a 4> a {\x,\ T t)d^ a {\y,\ T t) >~ A" 2 " 2 ^ < d a ^ x -i a (x,t)d a ^ x -i a (y,t) > . (40) 

The two point function can be evaluated in the continuum limit as A goes to infinity: 

< dMXx,\ T t)d a ij a {Xy,\ T t) >~ (Ax- Ay)" 2 " 2 ^. (41) 

Putting y = yields the large scale behaviour of the velocity two point function. The Fourier 
transform scales as 

E(k) ~ k 2d * +1 . (42) 

Notice that this is the infra-red spectrum. Once a fixed point is specified, the infrared spectrum 
follows. As expected, no prediction on the UV spectrum is derived. 

d) Some Examples 

Let us now give the simplest possible solution of the Navier-Stokes equation: the Gaussian 
free field with charges at infinity! 16 ]. Choose the Lagrangian to be: 

S = ^ j d 2 * W) 2 ( 43 ) 

and identify ip =: expz</> :. We explicitly suppose that there is a charge Q at infinity. This 
conformal field theory can be abstractly defined by its energy-momentum tensor and the operator 
product algebra. One can also use also a formal device to write down the action of the deformed 
Gaussian theory with a charge at infinity. Indeed, it is sufficient to add to the free Gaussian 
action a source term j d 2 xJ (X> (x)(J) where the source is = iQ5 2 (x — oo). This formal source 
at infinity implies that correlation functions of vertex operators are non-zero provided the sum of 
the exponents of each operator is equal to Q. Then the field ^2 =: exp2i0 : has dimension 2 + Q 
whereas = 0.5(1 + Q)- So vertex operators are good candidates to represent the asymptotic 
regime of the stream function. The charge at infinity is restricted to — 5 < Q < — 1. This theory is 
quite peculiar. Indeed, all UV properties, such as short distance expansions, are highly dependent 
on the value of the charge at infinity. In a sense, we can say that this reflects the non-unitarity 
of the theory and the increase of correlation functions with distance, i.e. large distance properties 
influence small scale characteristics. 

We would like to interpret this theory in a probabilistic way where operators are replaced by 
fields. This require that we identify : expicj) : with its formal adjoint : exp(— (1 + Q)i(f>) :. Therefore, 
we deduce that 2 : 

Q = -2. (44) 



The central charge is c = —11 
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Notice that in this case the infrared spectrum behaves as 

E ~ k°, (45) 

fe^O 

and in fact it happens this apparently innocuous solution is also turbulent. This demonstrates the 
idea that two dimensional turbulence is a problem where large scales are relevant. 



IV- Turbulent Fixed Points 



a) Currents in the Infra-Red Limit 

The conventional picture of turbulence requires that some physical entity (eg the energy) is 
transferred from large to small scales. This means that currents are non zero constants. We have 
proven that in the continuum limit n-enstrophies and the energy are conserved quantities. The 
conserved quantities can be written down as reciprocal-space integrals: 

H n , a = J dq h n , a {q) (46) 

The right-hand side is only modulus dependent. The integrand represents the conserved quantity 

per unit mode. As H n a is asymptotically conserved, there is a conserved current as well. More 
precisely, 

^ + d q R n , a (q) = (47) 
where R n ^ a is the current. The current involves the Fourier transform of u)%(x) and reads 

h nta {q) = 27tqL{- q ){u n a )(q) 

R nAv) = 4 / K,a( k ) (48) 

at Jk>q 

where L(q) is the Fourier transform of the characteristic function of a bounded domain of integra- 
tion L. The next step is to come back to real space. Introduce the structure factor 

9 q (x) = [ d 2 k e i2nk - x . (49) 



Then 

R d 

,11 . ^ 



" dt 



j d 2 x q *LU%- (50) 



The infinite volume limit can now be taken. Notice that the currents can be written as integrals 
of densities: 

Rn,a(<l) = / d 2 X r n , a (x) 

d (51) 

r n , a {x) = 6 q * -w^(k) 
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The time-derivative is explicitly taken as in the Euler equation. Currents appear as integrals of 
current densities which depend on q. In the Kolmogorov theory, current densities are constant and 
finite. We would like to evaluate current densities in the asymptotic regime. So rescale space and 
time such that 

r n>a (Xx, X T t, q) = A" d r n>A -i a (x, t) (52) 

where d is the dimension of r. The current density is given by a convolution; it can be simplified 
by noticing the scaling property: 

q (x-y) = X 2 6 j _(Xx-Xy). (53) 

Then 

r n ,\-i a (x,t,q) = -n j d 2 yX 2 0«.(Xx - Xy)e a/3 d a ip^a(y)di3^\-^a(y)^l-i a (y)- (54) 

We can now use the fact that, for large A, we have convergence of the approximate (^-function 
Qsl ~ 5 to obtain 

A 

r n ,\-i a (x,t) = -ne Q/3 a a ^ A -i a (x)a /3 AV'A-ia(^)^A r i 1 a( a; )- ( 55 ) 

This is ^-independent. So we have proven that in the asymptotic regime currents are q independent 
and represent local sources of turbulence. 



b) The Polyakov Equation 

Current densities are composite fields. Their field theoretic dimensions are the sum of two 
terms: the canonical dimension which is 2 as they are densities, and the anomalous dimension 
which accounts for the singularities appearing in short distance expansions. More precisely, the 
dimension of densities is given by 

d = 2 - d^ 2 - d^ n _^ + (n + l)d^. (56) 

The anomalous dimension cancels all divergences appearing in the short distance expansions of 

—j*- and cj^ _1 . Writing R n a (X T t) = j d 2 x A 2_d r n A -i a (x, t), we are interested in the asymptotic 

behaviour of the rescaled current densities r nt \-i a (x, t) = A 2_d r n A -i Q (x, t). The rescaled current 
density is a product of two fields: 

f BfA -i(M) = a d *-i + ^- 2 ^- 1 >-< n+1 >« l *A 2 t B - 1 )£^,A-i.V'n-i,A-i. (57) 

where L = (L_ 2 L 2 _ 1 - L_ 2 L 2 _ 1 ). The first power of A is the remnant of the derivatives which 
do not contribute to the anomalous dimension. This product is still divergent in the continuum 
limit. Only a fine tuning of the field dimensions can preserve a finite limit. Now, suppose that 
the leading part of short distance expansions at the fixed point defines an associative algebra, i.e. 
the leading part Dom(^) of field products satisfies Dom(^iDom(-02'03)) = Dom(Dom(-0iV'2)' i / ; 3) 
3 . This allows [V^HV'n-i] an d [V'HV'n] to be related, leading to 

r n ,\-i a (x,t) = A ^- 1 +^+ 2 ("- 1 )-^ + i fl rf Vi- 2 ("- 1 )-("+ 1 ) <i *^ Tl+u _ 1(l . (58) 

3 The operator algebra is always associative. If the stream function is a simple current^, then 
the leading part of short distance expansions defines an associative algebra. 
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The field on the right hand side is present as it is the lowest dimensional descendent field of y> ra +i 
which turns out to be a scalar (recall that Vn+i is pseudo-scalar). The right hand side converges 
provided 

d^ n+1 = d i)n _ 1 +d^ 2 +2(n- 1). (59) 

A solution can be easily found: 

<h n = \n(n + Q) (60) 

Notice that n = 2 gives the Polyakov equation. This generalised Polyakov equation was already 
mentioned in ref.fll]. We have proven that the current densities converge in the infrared limit. 
The scaling properties of the fields Cip n+ i imply that currents vanish: 

Rn = (61) 

This implies that the infrared behaviour of two-dimensional turbulence is not represented by a 
cascade. 

As a matter of fact, the Gaussian field theory with a charge Q = —2 at infinity satisfies (60). 
Moreover, vertex operators are simple currents; therefore the leading part of the operator algebra is 
associative. We have found a fixed point such that fluxes of n-enstrophies are finite in the infra-red 
limit. If the leading part of the operator algebra of the fixed point defines an associative algebra, 
the deformed Gaussian theory is the unique asymptotic description of 2d turbulence. 

V- The Coulomb Gas Picture 

a) The effective action 

The Gaussian fixed point fulfils all requirements to represent an asymptotic state of two- 
dimensional turbulence. Nevertheless, we would like to describe not only the fixed point but the 
effective theory in its basin of attraction. Recall that this theory is a perturbed conformal theory 
whose mass is getting smaller and smaller as one gets closer to the fixed point. 

Let us briefly sum up the time evolution of solutions of the Navier-stokes equation. First of 
all, there is a transient period when t < r where solutions evolve from the initial condition to a 
field theory characterised by an action S T . Two possibilities are to be envisaged. Suppose that 
a solution never falls within the basin of attraction of a conformal fixed point, then it can keep 
wandering in function space never reaching any steady state. This is an unattractive case. On 
the other hand, suppose that a solution falls within the basin of attraction of a fixed point after a 
transient period. We will show that the most general effective field theory in the basin of attraction 
of the Gaussian fixed point with a charge at infinity admits vortices. 

Let us now describe the effective field theory in the vicinity of the fixed point. This effective 
theory should respect all the symmetries of the fixed point. We assume that S T is in the basin of 
attraction of the deformed Gaussian fixed point. Let us notice that the deformed Gaussian fixed 
point is invariant under two discrete symmetries which leave the stream function invariant: 

4> — > —<j>, Q — » —Q 

(f> -> (f> + 2vr 

The first invariance ensures that is a real field whereas the second leaves the vertex operator 
invariant. The interaction term of the effective theory being an even periodic function can be 
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(62) 



expanded in a Fourier series: 



\ f 1 ~ > 

Seff = — I d 2 x{-d4> 2 a ~ Zn cosn M ( 63 ) 

77 J " n>l 

where ip a is now an effective field defined with a cut-off a. The interaction term can be separated 
in two terms, the first one containing relevant and marginal fields, the second containing irrelevant 
fields. In the infra-red regime, the deformed Gaussian theory is a renormalisation group fixed 
point. It can only be perturbed by irrelevant fields. The coupling constant of each irrelevant fields 
behaves as z n = zo n (X~ 1 a) 3 2 . Thus we end up with an effective potential: 

n(n-\-Q ) q — 

z Qn {a\- 1 )— 2 cosn</> a (64) 

n>n 

where n is the minimal integer such that n (n + Q) > 4; i.e. n = 4. In the infrared regime the 
effective potential goes to zero. The constants z§ n explicitly depend on the initial probability law. 
To understand the physical meaning of this effective theory, it is convenient to draw an analogy 
with a Coulomb gas. 

b) Turbulence as a Coulomb Gas 

We will show that the effective theory is equivalent to a neutral Coulomb Gas^ 16 ]. This is a 
gas of charged particles allowed to appear and disappear in a grand canonical fashion. This result 
stems from the expansion of the partition function: 

Z = Jv<j> eMSeff), (65) 

which can be written: 

i>n , mi,ni~0 1>6 i i<] 

where the product is not zero provided that the sum of all qi is —Q. This is nothing but the grand 
canonical partition function of a Coulomb gas where charges qi = — oo... — no, n ...oo interact 
according to a Coulomb potential: 

i>4 mi,ni=9i>no i i<j 

where V(x) = — ^ In |x| is the Coulomb potential. The gas is charged as the sum of the charges 
adds up to — Q. Each choice of m, particles of charge — qt, and of charge qi gives a factorial 
factor to take into account the undistinguishability of charges. The temperature of this gas is 
and particles of different charges have different fugacities Zj. Notice that the fugacities are very 
small. 

Let us now interpret the role of these charges for turbulence. The correlation functions of vertex 
operators are easily computed in the Coulomb Gas theory. Indeed, expz^ a (x) is represented by 
an insertion of a positive charge at point x. Correlation functions of vertex operators are given by 
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correlation functions of charges. Therefore, two-dimensional turbulence is equivalent to a Coulomb 
gas where the probabilistic properties of the stream function are given by the statistical behaviour 
of a positive charge. In view of numerical simulations, our aim is to determine the classical 
configurations of the stream function and the vorticity. To do so, we need a further modification 
of the Coulomb gas picture where the role of charges as vortices will be clear. 



c) Two Dimensional Turbulence as a Gas of Vortices 



We have just have shown that a charge q = 1 represents the stream function in a probabilistic 
sense, i.e. correlation functions of the stream function are equal to those of this charge in a Coulomb 
gas. Nevertheless, this does not give the final description of two-dimensional turbulence as a gas 
of charges. Indeed, we are about to see that these charges generate classical field configurations 
whose statistical behaviour is exactly the one expected for a Coulomb gas. 

Consider the 0(2) sigma model for stream function configurations, i.e. an action given by 

Seff = ^ J d 2 x{d^f (68) 

where the stream function ijj is now taken to be periodic. Notice that this is nothing but the kinetic 
energy of a configuration. As expected, the effective action can be expressed as a local action in 
the velocity field. The stream field is probabilistically almost a free field. Each stream function 
configuration is weighted with a Gaussian action. We allow singular stream functions: different sin- 
gularities correspond to sectors in the theory, and the total partition function is obtained summing 
over all sectors. The field -0 is a sum of a random part and a classical part: 

ip = 4> + A- (69) 
More precisely, the classical field is solution of the equation 

AVc = (70) 

i 

where defects of charge % have been considered and <5-functions are smeared over a distance a. A 
charge Q is put at infinity, the gas is neutral if all the charges add up to —Q. The classical solution 
reads 

^c(x) = y^&arctan 'f X ^ 2 . (71) 

. [X Xi)i 

The integer is the charge of the vortex. The random part <p will not contribute to defect 
correlation functions. This decomposition tells us that a stream function configuration is entirely 
specified by a countable number of vortices. 

The energy of each configuration is given by: 

S = ^- f d 2 x{d(t) + d^ c ) 2 . (72) 
4vr J 

This can be easily computed and gives: 

S=^([dMdcf>) 2 -^J2m^\^ l \)- (73) 

77 i<3 a 
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The energy is finite if and only if: 



J> + Q = o 



(74) 



This is nothing but the energy of a charged Coulomb gas. To complete this correspondence, let 
us choose to restrict charges to — oo... — n ,n ...oo. As already shown, the stream function is 
represented in the Coulomb gas picture by the insertion of a charge in the medium. It is even 
simpler in the 0{2) sigma model. In that case, the charge q x at x interacts with other free 
charges. The averaged response of the medium consists of a classical field ip c which depends on 
the fluctuating number of free charges. The correlation function of two charges representing the 
correlation function of the stream function is simply given by the average over the fluctuating 
number of particles of the weight exp — S(ip c (x, y)) depending on the position of the two fixed 
charges. This is nothing but the weight of the classical field i/) c (x,y). This result can be extended 
to an arbitrary correlation function . We conclude that these classical configurations represent the 
stream function of a two-dimensional turbulent flow. 

Finally, we can interpret the charge qi . Recall that the vorticity is minus the Laplacian of the 
stream field so that for the classical part we have: 



i.e. each configuration of two-dimensional turbulence is specified by a finite number of vortices, 
as shown by numerical simulations. The typical size of each vortex is a. Moreover, these vortices 
are quantised. Notice that the kinetic energy of a configuration is —2^^^ . qiqj In \xi — Xj\ and 
the n-enstrophies are proportional to ^ n the infrared regime, the probability of charge 

creation is zero so the energy as well as the n-enstrophies are zero. The convergence of energy and 
n-enstrophies to zero is with probability one. 

We have now proven that in the vicinity of the deformed Gaussian fixed point, the solutions of 
the Navier-Stokes equations are given by a gas of singularities where the stream function is a sum 
of vortices. The non-trivial infrared limit is due to the presence of an integer vortex at infinity. 
This provides a solvable example of two-dimensional turbulence. Moreover, the role of boundary 
conditions is made explicit. 



CONCLUSION 



We have shown that the asymptotic regime of two dimensional turbulence in the basin of at- 
traction of the deformed Gaussian fixed point is described by a Coulomb gas. This is an interacting 
neutral gas filled with charges \q\ > 4. We have explicitly shown that if the leading part of the 
operator algebra of the fixed point is associative, the deformed Gaussian fixed point is the unique 
asymptotic description of 2d turbulence. If universality arguments can be applied to self-organised 
criticality, it is clear that this must be the unique solution. 

Let us comment on the comparison with numerical simulations. First of all, it is conspicuous 
that the Coulomb gas has a natural description in terms of vortices. This is indeed what simulations 
obtain, i.e. coherent and stable vortices surrounded by a motionless flow. As initial conditions are 
always chosen to be Gaussian, we expect that as soon as defects appear, solutions fall within the 
basin of attraction of the Gaussian fixed point with charges at infinity. These defects are likely to 
be remnants from large initial fluctuations of the vorticity field. It is then natural to see solutions 




(75) 
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flow rapidly along renormalisation trajectories towards a Coulomb gas behaviour. To support this 
explanation, it would be relevant to know if the infrared energy spectrum is universal. We predict 
that the infrared energy spectrum is flat, i.e. behaves as k°. It would also be relevant to describe 
the early development of solutions and to prove that defects appear in a finite time. This would 
support our claim that the subsequent evolution of solutions is given by a Coulomb gas. 
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